Accuracy of the diffusion equation with extrapolated-boundary condition for 
transmittance of light through a turbid medium 
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The linear intensity profile of multiply scattered light in a slab geometry extrapolates to zero at 
a certain distance beyond the boundary. The diffusion equation with this "extrapolated boundary 
condition" has been used in the literature to obtain analytical formulas for the transmittance of 
light through the slab as a function of angle of incidence and refractive index. The accuracy of 
these formulas is determined by comparison with a numerical solution of the Boltzmann equation 
for radiative transfer. 
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I. INTRODUCTION 

Multiple scattering of light in a turbid medium is well 
described by the theory of radiative transfer [l|, 0, 3- 
This theory is based on a Boltzmann equation for the sta- 
tionary intensity I(f, s) of monochromatic light at posi- 
tion f and with wave vector in the direction s. In the sim- 
ple case of isotropic and non-absorbing scatterers (with 
mean free path I), the Boltzmann equation takes the form 



I s ■ Vl(r, s) = -I(r, s) + 1(f). 



(1) 



Far from boundaries the angle-averaged intensity, given 
here for 3 dimensions, 



1(f) = 



ds 
-in 



I(f, s) 



satisfies the diffusion equation 

V 2 /(r) = 0, 



(2) 



(3) 



which is easier to solve than the Boltzmann equation. 
The diffusion equation breaks down within a few mean 
free paths from the boundary, and one needs to return to 
the Boltzmann equation in order to determine 1(f) near 
the boundaries. 

A great deal of work has been done on the choice 
of boundary conditions for the diffusion equation which 
effectively incorporate the non-diffusive boundary layer 
0, H S IE IS 0, II • These studies have led to the so called 
"extrapolated-boundary condition" 



I(r) = -£n-VI(f), 



(4) 



where r is a point on the boundary and n is a unit vector 
perpendicular to the boundary and pointing outwards. 
Equation Q implies that a linear density profile extrap- 
olates to zero at a distance £ beyond the boundary. The 
extrapolation length £ is of the order of the mean free 
path. 

In this paper we consider transmission through a slab 
of finite thickness L. We compute the transmittance 



T, the ratio between the incident flux and the trans- 
mitted flux, by solving the Boltzmann equation numer- 
ically. Previous work on this problem used the diffu- 
sion equation with the extrapolated-boundary condition 
[E fiol fill. Il3 | (or an alternative large-L// approxima- 
tion |7j) to derive convenient analytical formulas for the 
dependence of T on L and I, on the refractive index of 
the slab, and on the angle of incidence. It is the pur- 
pose of the present study to determine the accuracy of 
these formulas, by comparison with the results from the 
Boltzmann equation. We generalize and extend a pre- 
vious study by De Jong [13] in the context of electrical 
conduction through a disordered metal, where the issue 
of refractive-index mismatch and angular resolution has 
not been considered. 



II. CALCULATION OF THE TRANSMITTANCE 
FROM THE BOLTZMANN EQUATION 

We consider a slab containing a disordered medium 
between the planes 2 = and z = L (see Fig. [IJ . The 
intensity I(f, s) depends only on the z-coordinate and on 
the angle between s and the z-axis. We define fj, = cos 0. 
The Boltzmann equation |T|) takes the form 



I(z) 



-I(z,n)+I(z), 



d/j,I(z, jj). 



< z < L,(5a) 
(5b) 



We supplement Eq. ([5]) with boundary conditions at Z = 
and z = L that describe reflection due to a refractive 
index mismatch, with reflection probability R(fJ-) : 



1(0, M )= R(ji)I(0,-fi) 
I(L,-iT)= R(n)I(L,n). 



Io((J>), 



H > 0, (6a) 
fj, > 0, (6b) 



The boundary condition at z — also contains the in- 
tensity due to a planar source with angular distribution 





-^ z 







FIG. 1: A sketch of the slab geometry. The disordered 
medium (thickness L) has a refractive index n (relative to 
the outside) and a mean free path I. The slab is illuminated 
by a plane wave incident at an angle So, vac, which is refracted 
to an angle #o inside the medium. The transmittance T is the 
ratio of the transmitted flux F and the incident flux Fq. 



7o(a0 inside the medium. Note that the angular distri- 
bution is different from that outside the medium, due to 
refraction at the boundary and due to the fact that part 
of the light is reflected before even entering the medium. 
The Boltzmann equation ([5]) with boundary condi- 
tions ([6]) implies for I(z) an integral equation of the 
Schwarzschild-Milne type [J, H, 0, [l3l] 



I{z)= M (z)+ dz' I(z') M 1 (z-z') + 
Jo L 

M 2 (z + z') + M 2 [(L-z) + {L~z')\ + 
M 3 (z - z') + M 3 [(L-z) - (L-z')] 

We have defined the kernels 

rf M c -l»l/te 
2Z// 

2l V 

1 ^ N(^R 2 ( f i)e-^ L+z ^ 1 ^. 



Mr(z) = 

M 2 (z) = 

M 3 (z) - 

/o 2i^ 

The factor N is given by 



NQi)= l-i? 2 (M)e 



-2L/l[i 



(7) 

(8a) 
(8b) 
(8c) 

(9) 



The kernels Mi , M 2 , and M3 describe propagation from 
z' to z with zero, an odd number, and an even number 
of reflections, respectively. The source term M is given 
by 



M {z) = \j dnN(»)I (n) 



,-z/lix_ 



■R{n)e 



-(2L-z)/lfi 



Once I{z) is known, the intensities I(z,ii) and I(z,—/j,) 
with ji > follow from 



I(z,u) = Io(»)N(n)e- z ^ + 



dz' 



(z'-z)/lii 



I(z') 



JV(M)%)e- ,/ "' 
L dz' ( 

JJI\ e 



-z /Ifl 



R(n)e 



< 2L - z 'V l »)l(z'), 
(11a) 



I(z,-fi) = I (n)N(n)R(n)e-V L - z V l <- 
^e^- z '^I(z') 



N{n)R(ji)e- ( - 2L -* )/l " 

i?(/i)e- z ' / ^)/(z'). (lib) 






This is a solution of the Boltzmann equation ^ with 
boundary conditions ([6]), as can be checked by substitu- 
tion. Integration over all /x then yields the Schwarzschild- 
Milne equation (J7J). We solve the integral equation Q 
numerically, by discretizing the interval (0, L), so that it 
reduces to a matrix equation [13j |. 

The quantity of interest is the transmittance T, defined 
as the ratio of the flux F that is transmitted through the 
slab and the flux Fq incident from the source, 



T = F/F . 
The transmitted flux is given by 



F = 2tt 



d(tnl(z,n), 



(12) 



(13) 



where c is the speed of light in the medium. The flux 
is independent of z, because there is no absorption. The 
total incident flux (including the flux which is reflected at 
the slab boundary before entering the medium) is given 
by 



Fn = 2n 






(14) 



We assume that the medium in the slab has a refractive 
index n > 1 (relative to the refractive index outside the 
slab). The lower bound \x c in the integral, defined by 
\x c = (1 — 1/n 2 ) 1 / 2 , is the cosine of the angle at which 
total internal reflection occurs [i?(/x) = 1 for /i < /x c ]. 
For < fi < fi c the reflection probability is given by the 
Fresnel formula for unpolarized light, 



Rfa) 



i-h 



//, v 



nfj, 



1 
Mvac + nn 2 
n 2 (l- M 2 )] 1/2 . 



nfMi 



/' 



n/M/ 



(10) 



The relation between /z vac = cos 9 vac and /x = 
Snell's law, such that the angle of incidence 9 V 



, (15a) 
(15b) 

cos 8 is 
outside 




0,vac 



FIG. 



2: Transmittance as a function of the incident angle 
#o,vac for two different values of L and n. The thick curves 
are computed from the Boltzmann equation, the thin curves 
are the diffusion approximation (I20|l , 



the medium (in "vacuum") is refracted to an angle in 
the medium (see Fig.QJ. 

We have calculated the transmittance for the case of 
plane- wave illumination, Iq (/x) — IoS(fi-fio)- The results 
are shown in Fig. [2 where T is plotted as a function of 
the angle of incidence #o,vac outside the medium [/Uo,vac — 
cos #o, vac is related to /x by Eq. (|15b[) ]. We show results 
for two different ratios L/l and two values of n (thick 
curves). For n = 1 the transmittance is non-zero for 
all incident angles, but numerical difficulties prevent us 
from going beyond #o. V ac — 87°. The thin curves in Fig. [2] 
are the results of the diffusion approximation, which we 
discuss in the following section. 



III. COMPARISON WITH THE DIFFUSION 
APPROXIMATION 

The diffusion approximation for the transmittance has 
been studied by several authors [3, [l(J, [CO, Q3 • Here we 
briefly describe this approach, and then compare the re- 
sult with the numerical solution of the Boltzmann equa- 
tion. 

In a slab geometry the diffusion equation with 
extrapolated-boundary condition takes the form [cf. 
Eqs. © and ©] 



with boundary conditions 

j(o) = e/'(o), 



< z < L, (16) 



I(L) = -£/'(L). (17) 



We assume plane- wave illumination of the boundary z = 
0, at an angle #o.vac with the positive z-axis. A fraction 



1 — R(no) of the incident flux Fo enters the medium, 
and is first scattered on average at z = /x L [We recall 
that jUo = cos ^o, where 8 corresponds to #o,vac after 
refraction, cf. Eq. (|15b[) .] This plane- wave illumination is 
incorporated into the diffusion equation (| 16[) as a source 
term, 

i /(z) + id [l ~ R (^ Fo6 ( z - W>0 = 0. (18) 
The solution of Eq. (fl8|) with boundary condition (JTTJ) is 



I(z) 



{ 3 te + z)(L + g- w Q 
4ttZ L + 2£ 

x [1 - R(fjo)]F , if < z < fjtol, 

3 (g + ^ Q(£ + e-z) 
4tt/ L + 2^ 

x [1 - R(fio)]Fo, if jj, l < z < L. 



(19) 
The transmitted flux F = —jTrl'(L) divided by the inci- 
dent flux Fq leads to the transmittance in the diffusion 
approximation, 



T diff = [1 - R(no)} 



£, + vol 



(20) 



This simple analytical formula combines results in the 
literature by Kaplan et al. [9] (who considered normal 
incidence) and by Nieuwenhuizen and Luck [7j (who con- 
sidered the Schwarzschild-Milne equation in the diffusive 
limit L > I). 

We still need to specify the value of the extrapolation 
length £. We will use an expression due to Zhu, Pine, 
and Weitz [11] , 



£ 



2 1 + C a 

3 1-d' 



(21) 



where the coefficients C\ and Ci are the first two mo- 
ments of i?(/x), 

d = 2 / dnnR{n), (22a) 

Jo 

C 2 = 3 / dnn 2 R{^)i (22b) 

Jo 

normalized such that C± = C2 = R for an angle- 
independent reflection probability R(n) = R. Compar- 
ison of Eq. (|2"Tj) with a numerical solution of the Boltz- 
mann equation in a semi-infinite medium by Aronson 8] 
shows that it accurately describes the length over which 
the linear density profile extrapolates to zero. The differ- 
ence is largest for n = 1, when Eq. ([2T|) gives £ = |Z while 
the Boltzmann equation gives an extrapolation length of 
0.7104/ which is somewhat larger [l|, 0|. The transmit- 
tance Tdiff is compared in Figs. [2] and [3] with the exact T 
from the Boltzmann equation. 

Once the transmittance T for plane-wave illumination 
as a function of /uo, V ac = cos #o. va c is known, one can 
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FIG. 3: Ratio of the transmittance Tdig according to the 
diffusion approximation and the exact result T according to 
the Boltzmann equation, for normal incidence So, vac — 9q = 0. 
The inset shows T as a function of L/l, for the same values 
of n as the main plot. 



alytical formula 



Tdiff . 



3L 

IT 



i + c 2 

l-Ci 



(24) 



In the absence of a refractive-index mismatch (n = 1, 
C\ = C-2. = 0) this formula has been found [l3| to agree 
with the Boltzmann equation within 3% for all L/l. For 
n > 1 the relative error in Eq. (|24|) is comparable to 
that shown in Fig. [3] for the transmittance at normal 
incidence. 

In conclusion, we have computed the transmittance of 
a turbid medium of mean free path I and length L from 
the Boltzmann equation as a function of the angle of 
incidence. We compared the results from the diffusion 
equation to this exact solution. The difference between 
the two transmittances stays below 6% for L > 31 and 
1 < n < 2. The diffusion approximation overestimates 
the transmittance for n = 1 and underestimates it in 
the presence of a significant refractive index mismatch. 
The relative error is largest for large refractive index mis- 
match. 



compute the transmittance Ttot for diffusive illumination 
by integrating over the angles of incidence, 

Ttot — 2 / d/iO,vac MO.vac ^(/io.vac)- (23) 

Jo 

The diffusion approximation (|2H)l and (|2~lj) yields the an- 
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